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We present an alternative for the description of galactic halos based on Tsallis' 
non-extensive entropy formalism; on this scheme, halos are stellar polytropes char- 
acterized by three parameters, the central density, p c , the central velocity dispersion, 
a c and the polytropic index, n. To evaluate these parameters we take the Navarro- 
Frenk- White paradigm as a comparative model and make the following assumptions: 
both halo models must have the same virial mass, the same total energy and the 
same maximal velocity. These three conditions fix all the parameters for a given 
stellar polytrope allowing us to compare both halo models. The halos studied have 
virial masses on the range 10 12 — 10 15 Mq, and it was found after the analysis that 
they are described, at all scales, by almost the same polytropic index, n ~ 4.8, im- 
plying an empirical estimation of Tsallis non-extensive parameter for this type of 
dynamical systems: q ~ 1.3. 

PACS numbers: 



INTRODUCTION 



Cold Dark Matter (CDM) models based on N-body numerical simulations predict exces- 
sive substructure and cuspy dark matter halo profiles that are not observed in the rotation 
curves of dwarf and LSB galaxies The significance of this discrepancy with observa- 
tions is still under dispute, leading to various theoretical alternatives, either within the 
thermal paradigm (self-interacting [2f and/or "warm" [3] DM, made of lighter particles), or 
non-thermal dark matter (DM) models (real 4] or complex [5] scalar fields, axions, etc). 
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However, none of these alternatives is free of controversy. On the other hand, the CDM 
model of collision-less WIMP's could remains still as a viable model to account for DM in 
galactic halos, provided there is a mechanism to explain the discrepancies of this model with 
observations in the center of galaxies. The main goal of this paper is to give an approach 
for such a possibility. The main idea for this alternative is as follows. Since gravity is 
a long-range interaction and virialized self-gravitating systems are characterized by non- 
extensive forms of entropy and energy, it is reasonable to expect that the final configurations 
of halo structure predicted by N-body simulations must be, somehow, related with states 
of relaxation associated with non-extensive formulations of Statistical Mechanics. 

The usual statistical mechanic treatment of self-gravitational systems is provided by the 
micro-canonical ensemble, which is compatible with negative heat capacities associated with 
known effects, such as gravothermal instability An alternative formalism that allows 

non-extensive forms for entrop y a nd energy has been developed by Tsallis |8( and applied to 
self-gravitating systems 0,[lO,lllj], under the assumption of a kinetic theory treatment and a 
mean field approximation. As opposed to the Maxwell-Boltzmann distribution that follows 
as the equilibrium state associated with the usual Boltzmann-Gibbs entropy functional, 
the Tsallis functional yields as equilibrium state the "stellar polytrope" , characterized by a 
polytropic index n. The stellar polytrope yields a Maxwell-Boltzmann distribution function 
(the isothermal sphere) in the limit n — > oo. This index is related to the "non-extensivity" 
parameter q of Tsallis entropy functional, so that the "extensivity" limit q — ► 1 corresponds 
to the isothermal sphere. 

Although the self-gravitating collision-less and virialized gas that makes up galactic halos 
is far from the state of gravothermal catastrophe, it is reasonable to assume that it is near 
some form of relaxation equilibrium characterized by non-extensive forms of entropy and 
energy. On the other hand, high precision N-body numerical simulations, perhaps the most 
powerful method available for understanding gravitational clustering, roughly yield density, 
mass and rotation velocity profiles that seem to fit observed galactic halo structures (pending 
the controversy on excess substructure and cuspy density cores specially on LSB galaxies). 
Admitting that stellar polytropes follow from an idealized approach based on kinetic theory 
and an isotropic distribution function, it is still interesting to verify empirically if the struc- 
tural parameters of the halo gas, or at least of the outcome of numerical simulations, can 
be adjusted to those of a stellar polytrope, the equilibrium state under Tsallis' formalism. 
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It is important to notice that the main objective of this paper is not to compare this 
polytropic model of dark halo with observational results coming from actual galaxies but 
with the Navarro-Frenk- White paradigm of dark halos which adequately describes the 
rotation curves of most galaxies. It is known that the NFW profile fits well the density profile 
of galaxies in the region outside their core. It fails in the central regions where observations 
show that the density profile is almost flat. In this work we show that, for example, the best 
polytropic fit to halos with NFW profiles follows from polytropes characterized by densities 
in the range 3.7 x lO~ 4 M /pc 3 < p c < 1.2 x 10 _3 M Q /pc 3 and polytropic indices is almost 
constant with a value near n = 4.8, or q m 1.3 

Since these best-fit polytropes have the same observable quantities as the NFW halos 
without central density cusps, they might provide an even better fit to halo structures than 
the usual NFW profiles. Furthermore, the present analysis and results can be used to 
calibrate the values of the free parameter q that emerges from Tsallis's formalism. 

Hence, we propose to verify which parameters of the stellar polytropes provide a suitable 
description of the halo that resembles the one that emerges from the well known numerical 
simulations of Navarro-Frenk- White. For the wide virial mass range of 10 10 < M vir /M < 
10 15 . 

The paper is organized as follows: in section II we provide the equilibrium equations 
of stellar polytropes and briefly summarize their connection to Tsallis' non-extensive en- 
tropy formalism. In section III we discuss the dynamical variables of NFW halos, while in 
section IV we describe a procedure to compare a polytropic halo with an NFW one and ob- 
tain numerically the parameters that characterize such polytropic halos, producing graphics 
showing such comparison. A summary of our results is given in section V. 



II. TSALLIS ENTROPY AND STELLAR POLYTROPES. 



For a face space given by (r, p_h the kinetic theory entropy functional associated with 



Tsallis' formalism is |9|,ll0|,and[ll| 



S q = ——[ J (/WK^p, (1) 

where / is the distribution function and q > 1 is a real number. In the limit q — > 1, the 
functional (0) leads to the usual Boltzmann-Gibbs functional. The condition 5 S q = leads 
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to the distribution function that corresponds to a stellar polytrope characterized by the 
equation of state 

p = K nP l+1 /\ (2) 

where K n is a function of the polytropic index n, and can be expressed in terms of the 
central parameters: 

Kn = -^7-. (3) 

p x/n y > 

The polytropic index, n, is related to the Tsallis' parameter q > 1 by: 

3 1 

n =2 + —l (4) 
The stability condition 5 2 S g < is only satisfied generically for polytropes with 3/2 < n < 5 
(or 9/7 < q), which are then stable equilibrium configurations. Polytropes with n > 5 are 
then meta-stable configuration which undergo gravothermal instability for sufficiently large 
density contrast p/p c (see jlO, 11]). 



The standard approach for studying spherically symmetric hydrostatic equilibrium in 
stellar polytropes follows from inserting (J2J) into Poisson's equation, leading to the well 
known Lane-Emden equation jl^] 

2^^ + *" = o, (5) 

with 



x 2 dx \ dx J 

. l/n 
P_ X 

Pc 



(6) 



r _ 2 AnGp c 2 p c 

x = v r ° = ^T' ac = (7) 

G = 4.297 x 10- 6 (km/ ^ c)2 (8) 
M /kpc 

where a c and p c are the central velocity dispersion and central mass density respectivelly; 
and we take this value for the gravitational constant due to the units we are using. Notice 
that the velocity dispersion is a measure of the kinetic temperature of the gas by the relation: 
a\ = k B T c /m, with k B being Boltzmann's constant, and that we are using a normalization 
for t~o, which differs from the usual one by a factor l/(n + 1). We find it more convenient to 
consider instead of (jHJ) the following set of equivalent equilibrium equations 



dx 

dZ n M Z 1 - 1 ^ 

dx n + 1 x 



2 • 0) 
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where M. and Z relate to M, p (mass and mass density at radius r) by 

M - =£r 2 -i (10) 

Notice that in the limit n — > oo, equations (J^J) become the equilibrium equations of the 
isothermal sphere. 

Once the system (jOJ) has been integrated numerically, the velocity profile derived from 
the virial theorem takes the form: 

V 2 (x) = (11) 

where a c can be given in km/sec. The radial distance r in kpc and enclosed mass M(r) in 
solar masses are given (from in terms of x and M. by 

r/kpc = 0.004220 -^- f Mo/pc^ 172 
km/sec \ p c J 

M/M = 944.9737 (V^— Y (^^Ivf)^ M , (1 2 ) 



km/secy \ p, 

Another important dynamical quantity is the total energy of the stellar polytrope 

E = K + U 

= - / 47rr 2 P r dr - dr- y ' y ' 



2 Jo Jo r dr 

(13) 



leading to 

-3 CM 2 /p\V™ 



poly 



n — 5 



2 r Vp, 



(-) n a c (-(n + l)M„-(n-2)4 7rr 3 Pl /)] (14) 



which must be evaluated at a fixed, but arbitrary, value of r marking a cut-off scale. 



III. NFW HALOS 



s yield the fol 
HQ, Q: 



NFW numerical simulations yield the following expression for the density profile of viri- 
alized galactic halo structures 



y (l + vY 



»0 HO r\ 
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where: 



3 [In (1 + c ) -c /(l + c )]' 



p crit n h 2 = 253.8 h 2 fi 



M, 



© 



kpc 



3 ' 



y = c —, 



(16) 

(17) 
(18) 



where f2 is the ratio of the total density to the critical density of the Universe, being one for 
a flat Universe. Notice that we are using a scale parameter (y) that is different from that of 
the stellar polytropes (x). The NFW virial radius, r vir , is given in terms of the virial mass, 
M vir , by the condition that average halo density equals A times the cosmo logical density p : 



Ap 



3M vir 



(19) 



4 7rr 3 

vir 

where A is a model-dependent numerical factor (for a ACDM model we have A ~ 100 
). It is important to remark that this last relation between the mass and 
virial radius in terms of cosmological parameters, Eq. (jl9j) . is valid for any halo model. The 
concentration parameter c can be expressed in terms of M vir by 1?1 13| : 



c ~ 62.1 x 



M Yir h 



-0.06 



(20) 



hence all quantities depend on a single free parameter M vir . The mass function and circular 

r2 f „\ _ 4 n GM/r, leading to: 
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velocity follow from (|15|) by M(r) = Arc J pr 2 dr and V 2 (r) 



M = 47T 



<Wo 



y 



V 2 =4nG5 oPo ^ 
Co 



Ml + 2/)" , 

1 + y 

Mi + y) 1 



(21) 
(22) 



y 1+2/. 

Given p and M, the gravitational potential per unit mass and radial and tangential pressures 
follow from: 



which combine to give: 



P ' 

r 

P ' 



G — , 



2{P r ~ Pj.) 



-G 



Mp 2a{r) 



(23) 
(24) 

(25) 
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FIG. 1: Mass profiles for a NFW halo with M vir = 10 12 M & and r vir = 260 kpc (solid curve) and 
compared fit stellar polytrope (dashed curve). 

where: 

a = (26) 

is the anisotropy parameter < |a| < 1, so that a = corresponds to an isotropic velocity 
distribution, since pressure is proportional to velocity dispersion. The parameter a is often 
taken as a constant in the range < a < 0.2, or given by the more realistic empirical ansatz 
of Ostipov and Merritt QQQ- 

The total energy for the NFW halo follows from the general expression (1131) . where P 

nn 

should be obtained from the integration of (|23|) for a specific form of a. As shown by |15l.l2l|. 
the curves of P = pa 2 obtained from the Ostipov-Merritt ansatz are very close to those 
of the isotropic case (a = 0), hence we will consider only isotropic NFW halos. Although 
analytic solutions of (}25|) exist for a = (see [lill^]]), we will use instead the evaluation of 



(|13|) obtained by Q| in which the leading term of the total energy is given by: 

Mj 

Enfw = — G v Fo (27) 
2r vir 

where F Q has the approximate values: 

2 , / c 



Ft 







3 V21 

and Cq is given in terms of the virial mass by (|2*U|) 



k) ' < 28 » 
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IV. POLYTROPIC AND NFW HALOS COMPARISON. 

In order to compare stellar polytropes to NFW halos, it is important to make various 
physically motivated assumptions. First, we want both models to describe a halo of the 
same size, but since the virial radius, r vir , is the natural "cut-off" length scale at which the 
halo can be treated as an isolated object in equilibrium, "same size" must mean same virial 
mass, M vir , by equation (fTTJ|h 

Secondly, both models must have the same maximal value for the rotation velocity ob- 
tained from (|TT| and (|22|) . This is a plausible assumption, as it is based on the TuHy-Fisher 
relation, j^J , a very well established result that has been tested successfully for galactic sys- 
tems, showing a strong correlation between the total luminosity of a galaxy and its maximal 
rotation velocity. It can be shown that the Tully-Fisher relation has a cosmological origin, 



|17| , associated with the primordial power spectrum of fluctuations (the so called "cosmolog- 
ical Tully-Fisher relation"), hence it is possible to translate the correlation between maximal 
rotation velocity and total luminosity to a correlation between maximal rotation velocity 
and total (i.e. virial) mass. Since, by construction we are assuming the polytropic and NFW 
halos to have the same M vir , their maximal rotation velocity must also coincide. 

Our third assumption is that the polytropic and NFW halos, complying with the previous 
requirements, also have the same total energy computed from ^Tj\ and from (fTlj) evaluated 
at the cut-off scale r = r vir . The main justification for this assumption follows from the fact 
that the total energy is a fixed quantity in the collapse and subsequent virialized equilibrium 
of dark matter halos j^j]. 

Since all structural variables of the NFW halo depend only on the virial mass, once 
we provide a specific value for M vir all variables become determined in terms of physical 
units by means of the scale equation ()18|). Polytropic halos, on the other hand, lack a closed 
analytic expression for mass, velocity and density profiles. In this case, equations (or (jo^)) 
yield numerical solutions for these profiles expressed in terms of the three free parameters 
{p c , a c , n}. The comparison of these profiles with those of the NFW halos requires that we 
find explicit values of these free parameters, so that the conditions that we have outlined 
are met. Since we have selected three comparison criteria for three parameters, we have a 
mathematically consistent problem. The first constraint on {p c , a c , n} follows by demanding 
that, for a given M vir characterizing a NFW halo, we have M(r vir ) = M vir for the polytropic 
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TABLE I: Parameters characterizing the polytropes while being compared to NFW halos 
log lo (M vir /M ) p c [M /pc 3 ] a c [Km/s] n q K n v mSLX [Km/s] r vir [kpc] 



15 


3.7 x 10" 


-l 


982 


4.93 1.29 4873.4 


1504 


2606.2 


12 


7.5 x 10" 


-4 


108 


4.87 1.30 478.94 


164 


260.6 


11 


9.0 x 10- 


-4 


52 


4.83 1.30 221.82 


79.1 


120.9 


10 


1.2 x 10- 


-3 


25 


4.82 1.30 100.68 


38.2 


56.1 



halo obtained from the numerical solution of (jHJ) with M and r given by (|12j) . A second 
constraint on the polytropic parameters follows from equating E po i y from (JHJ), evaluated at 
r = r vir and M = M vir , with -Enfw from ([27)1 . Having fixed two of the polytropic parameters, 
the third one can be fixed by demanding same maximal velocities in the curves for and 

(IP- 
Following the guidelines described above, we proceed to compare NFW and polytropic 
halos for M vir ranging from 10 10 up to 10 15 solar masses. From the present comparison we 
find that the values for central density, p c , of the polytropic halos are inversely proportional 
to M vir , while the values for the central velocity dispersion, a c , are directly proportional to it 
(this is expected, since a c is a scale parameter in self-gravitating systems). The polytropic 
index, n, is almost constant for the selected range of M vin showing a very small growth as 
M vir increases. This implies the same qualitative behavior of the Tsallis parameter q: it 
is also almost constant and is slowly increasing as the virial mass grows. It is worthwhile 
mentioning that the proportionally term K n in the polytropic equation of state (JBJ) shows 
a very noticeable change, rapidly growing as M vir increases. All these results are displayed 
explicitly in table |TJ The figures depict the mass profiles, figure [TJ the velocity profiles, 
figure El and the density profiles, figure El for the resulting polytropes with M vir = 10 12 M , 
juxtaposed with the same profiles for a NFW halo with same M vir . For other values of M vir 
the mass, velocity and density profiles are qualitatively similar to the ones displayed in these 
figures. 

V. CONCLUSIONS 

Motivated by the fact that stellar polytropes are the equilibrium state in Tsallis' non- 
extensive entropy formalism, we have found the structural parameters of those stellar poly- 
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FIG. 2: Velocity profiles for the same NFW halo of figure 1 (solid line) and its compared stellar 
polytrope (dashed curve). 
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FIG. 3: Density profiles for the same NFW halo of figures 1 and 2 (solid curve) and its compared 
stellar polytrope (dashed curve). 

tropes that allows us to compare them with NFW halos of virial masses in the range 
10 10 < M vir /M < 10 15 . The criteria for this comparison consists in demanding that the 
polytropes describe a halo having the same virial mass, virial radius, maximal rotation ve- 
locity and total energy as the NFW halo. These three conditions are sufficient to determine 
the three structural parameters {p c , cr c , n} of the polytropic model; the results are displayed 
in Table |U 
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It is important to emphasize that the criteria that determine these polytropes are based 
on physically motivated assumptions: the virial radius and mass are the natural parameters 
characterizing the size of a given halo, same maximal velocity follows from the Tully-Fisher 
relation, while same total energy follows from the virilization process. As shown in Figure 
the mass distribution of the polytrope grows much slower than that of the NFW halo 
up to a large radius (100 kpc) containing the core and the region where visible matter 
concentrates. Hence, as shown by Figure |2l the velocity profile of the polytrope is much 
less steep in the same region than that of the NFW halo. These features are consistent 
with the fact that NFW profiles predict more dark matter mass concentration than what 
is actually observed in a large sample of galaxies 0, 24,13]. Also, as shown in Figure 01 
the obtained polytropes have flat cores, very similar to the flat isothermal cores observed 
in LSB galaxies (as a contrast, the cuspy cores of NFW halos seem to be at odds with 
these observations (26, 27, 2sj], also [s^Ql)- This flat density core is a nice property, which 
combined with reasonable mass and velocity profiles, qualifies these polytropes as reasonable 
(albeit idealized) models of halo structures. 

However, in spite of their nice theoretical properties (i.e. their connection to Tsallis' for- 
malism) and reasonable similarity with equivalent NFW halos, the stellar polytropes we have 
examined are very idealized configurations and so we are not claiming that they provide a 
realistic description of halo structures. Instead, we suggest that their described features and 
their connection with Tsallis' formalism might indicate that the latter could yield useful in- 
formation in understanding the evolution and virialization process of dark matter. Although 
it is necessary to pursue this idea by means of more sophisticated methods, including the 
use of numerical simulations along the lines pioneered by JjJ , the simple approach we have 
presented has already given interesting results. For example, with respect to the parameter, 
q, we recall that it is a free parameter of the Tsallis' non-extensive thermodynamics and 
which has not been fixed for the cosmological case. In this work, by using such statistic in 
cosmological systems, we are able to determine its behavior as a function of the virial mass, 
and turns out to be almost constant, with a values around q ~ 1 .2. This resu lt could be 



Tsallis 2001 1 . 



used in other contexts where the extended statistic is also applied 

It is well known |l2j that stellar polytropes with n < 5 (like King halos) have a finite 
cut-off scale and finite total mass, though for the polytropes that we have studied this 
cut-off scale is much larger than r vir (just as the "tidal radius" of King halos is much 
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larger than their virial radius). However, as shown in jloL 11 1, polytropes characterized by 
this polytropic index correspond to stable equilibrium states that are generically free from 
undergoing gravothermal instability. 

As mentioned, the results presented in the present work show that a dark matter halo 
made out of matter which satisfies a polytropic like equation of state, describes the halo in 
a way as good as the description obtained from the NFW numerical simmulations, that is 
their paradigm. Furthermore, our description is even nicer as long as it does not have great 
density growths near the center. However, these results does not directly imply that the dark 
matter halo do obey a non-extensive entropy formalism. Further tests and expriments are 
needed in order to consider that such formalism is the one describing the thermodynamics 
of actual dark matter halos. At the moment, this idea is a possibility which is reinforced by 
our analysis. We believe these properties to be very encouraging and are currently engaged 
in a more detailed examination of these polytropes 2^| . 
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